Reply to the comment on "Incomplete equilibrium in long-range interacting systems 
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After the rejection of their comment arXiv:cond-mat/0609399vl] to our Phys. Rev. Lett. 97, 
100601 (2006), the Authors informed us that an extended version of their comment is going to be 
published in a different journal under the direct editorial responsibility of one of them. We then 
decided to make publicly available our formal reply, originally prepared for publication in Phys. 
Rev. Lett. 

PACS numbers: 05.20.-y, 05.70.Ln, 05.10.-a 
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In their comment [l[ , Tsallis et al. present three points 
which they claim to confute the conclusions of our work 
in Ref. [2|. Here we show that this is not the case. 

An important issue about the quasi-stationary states 
(QSS's) displayed by isolated long-range interacting sys- 
tems before relaxing to equilibrium is whether the QSS's 
survive to the perturbations introduced by a thermal 
reservoir (canonical QSS's). In order to answer this ques- 
tion, we proposed [3] an Hamiltonian setup in which the 
microcanonical conditions are recovered when the cou- 
pling constant between system and thermal bath van- 
ishes. Indeed, we showed Q that canonical QSS's exist, 
with a lifetime which decreases as the coupling strength 
increases. In [2| we discuss the statistics of the system 
energy fluctuations and, contrarily to what is claimed 
in we do not ". . . extrapolate the conclusions for the 
canonical QSS's ...to the microcanonical ones", where 
by definition the energy fluctuations are zero. Hence, the 
first criticism raised in [l| offers to the unaware Reader 
a false representation of our results and motivations. 

On the other hand, we take here the opportunity of 
pointing out that the dynamical behavior during both mi- 
crocanonical and canonical QSS's is such that the inter- 
particle correlation is negligible (see Fig. [I}. Thus, in 
both cases the statistical mechanics description of the 
QSS's has to be based on the assumption of independence 
among elementary components and on the consequent 
application of the central limit theorem. This is consis- 
tent with our approach in where we give compelling 
evidence of the applicability of Kinchin's derivation of 
statistical mechanics 3 , once an appropriate estimation 
of the density of states uj(E) for the QSS's is given. It is 
also pertinent to add that the evidences we are present- 
ing in Fig. [1] provide full justification to the use of the 
Vlasov theory for the QSS's, as it is done, e,g., in 

The statement [l| that the non-Gaussianity of the 
one-body angular momentum PDF "... excludes the 
Boltzmann-Gibbs exponential form as the energy distri- 
bution in full phase space T . . . " is wrong. The Authors 
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FIG. 1: Particle-particle correlation during a microcanonical 
simulation of the Hamiltonian system in [3] with M = 5000 
particles. C x (x = 0,1) is defined as C x = ({xiXi+i) — 

~ (*i> 2 )«*? + /> ~ <^+/) 2 ) (Ml line) and 

C x = {xiXi+i) I '{^j (Xi)(x? +I )) (dashed line), with (■) = 

(X/i=i ")/! an d I < M/2. Since at t = all angles 9i's coin- 
cide (initial magnetization is 1), in (a) the dashed line starts 
from eg = 1. In (b) the angular momenta k's are initially uni- 
formly distributed and the dashed and the full lines coincide 
for all times. In all simulations we observed that C x during 
the QSS is not appreciably different than at equilibrium. The 
same is valid also for correlations of order higher than 2, even 
in canonical simulations or if the correlations are defined by 
averaging over different dynamical realizations. 



of [l[ are trivially proving that the PDF of a microstate 
(h,@i)i=i,2,...,M is not proportional to H " i,e ^ , H be- 
ing the Hamiltonian of the system with M particles. 
However, this does not imply that the total energy PDF 
must be different from p(E) — Ld(E)e~ l3E /Z, where (3 is 
an inverse temperature. On the contrary, in force of the 
central limit theorem p(E) still has the usual Boltzmann- 
Gibbs form Q , as we clearly verify in Q ■ The nontrivial 
point is that ui(E) for the canonical QSS's is a nonequi- 
librium density of state which can be calculated by con- 
sidering a submanifold of the T-space at constant mag- 
netization 0. 

In their final remark, the Authors of [l[ propose a non- 
exponential fitting of our results for \n[p(E)/oj(E)] after 
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the introduction of an ad hoc energy-shift E i— > E — 692 
for which they provide no explanation. A plot of the non- 
exponential function in [l[ without this energy-shift does 
not agree at all with the data from the dynamical simu- 
lations. We also point out that by fitting ln[p(E)/oj(E)] 
one implicitly assumes the validity of our calculation of 
uj(E), which is based on the fundamental thermodynamic 



relation linking temperature to the Boltzmann expres- 
sion for the entropy, S = kg la[u)(E)]. As we explain in 
Q, to propose an alternative to the exponential weight 
e~@ E without appropriately changing the determination 
of oj(E) is logically inconsistent. 

In summary, we have shown that the statements in [l| 
either are wrong or do not apply to our results in Q ■ 
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